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1.1  DESCRIPTION  OF  THE  PROBLEM 


The  results  of  an  experimental  and  theoretical  investigation  of  the  | 

I 

recirculating  flow  fields  induced  by  ducted,  variable-density,  turbulent  j 

jet  mixing  are  presented.  Recirculating  flow  fields  established  by  j 

turbulent  jet  mixing  of  two  coaxial  streams  in  a constant-area,  axisym-  j 

metric  duct  occur  in  many  Industrial  and  aerospace  burner,  furnace,  j 

or  combustor  configurations.  The  so-called  "sudden-expansion"  or 
"dump"  combustors  used  in  ramjet-rocket  propulsion  systems  are  designed  ' 

on  the  principle  of  establishing  and  maintaining  combustion  in  regions 
of  recirculating  flow  within  the  combustors. 

Figure  1 shows  the  essential  features  of  the  type  of  recirculating  j 

flow  field  studied  herein.  For  a certain  range  of  fluid  influx  condi- 
tions, ducted  jet  mixing  of  coaxial  streams  leads  to  the  creation  of 
an  eddy  of  recirculating  fluid,  existing  on  a time-averaged  basis  as 
a toroidal,  highly  vortical  region  with  high  levels  of  turbulent  in-  ^ 

tensities  and  relatively  low  average  velocities.  The  eddy  of  the 
recirculating  region  is  generated  when  the  "entrainment  capacity"  of 
the  higher  velocity  primary  jet  stream  is  not  satisfied  by  the 
weight  flow  of  the  outer,  lower  velocity  secondary  stream.  For 
those  conditions,  the  primary  jet  "establishes"  the  recirculation 
eddy  which  serves  to  feed  mass  into  the  jet  thus  dissipating  the 
jet  mechanical  energy. 

In  all  essential  time-averaged  features,  this  recirculating  flow 
is  a stationary,  separated,  turbulent,  toroidally  shaped  region  with 
and  eddy-like  structure  wherein  back-flowing  currents  carry  mass, 
energy,  and  momentum  from  downstream  toward  the  inlet  station.  The 
recirculation  eddy  can  exist  with  or  without  secondary  flow.  How- 
ever, at  a secondary  flow  rate  equal  to  or  greater  than  a critical 
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"blow-off"  value,  the  recirculation  eddy  structure  cannot  form,  and 
the  duct  flow  will  be  characterized  everywhere  by  positive  axial 
velocity  components.  The  pertinent  features  of  the  flow  field  with 
recirculation  are: 


1 . The  axial  location  of  the  two  points  on  the  duct 
wall,  denoted  and 
characteristic  time-averaged  stagnation  points  of 
the  flow, 

2.  The  time-averaged  axial  velocity  field, 

3.  The  tJme-averaged  distribution  of  mass,  energy 

or  temperature  of  the  fluid  from  the  primary  flow 
which  characterizes  the  qualitative  and  quantitative 
nature  of  the  mixing  that  has  occurred  between  the  two 
streams,  and 

4.  The  time-averaged  axial  distribution  of  wall  static 
pressure. 


in  Fig.  1,  which  are 


All  of  these  gross  features  of  the  flow,  their  spatial  distribu- 
tions and  magnitudes,  appear  from  experimental  evidence  to  be  depend- 
ent upon  the  ratios  (ug  /u  , P^/p  • )»  the  Reynolds  number 

based  on  the  primary  flow  at  the  nozzle  exit  station  (Nj^g  ),  and  upon 

p 

chemical  reactions,  if  any,  occurring  in  the  flow.  However,  it 
is  also  an  empirical  observation  that,  for  values  of  the  primary 
jet  Reynolds  number  defined  by 


N 


p u R 
P P P 


(1) 
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greater  than  about  1.5  x 10  , the  time-averaged  velocity,  species 
concentration,  and  static-pressure  fields  become  Independent  of  the 
jet  Reynolds  number  (%ep) » this  condition  is  satisfied  in  the 
present  study.  In  the  present  investigation,  the  flow  fields  studied 
were  nonreactive,  hence  of  uniform  total  temperature.  However, 
density  variations  were  created  throughout  the  flow  fields  by  the 
mixing  of  jets  of  different  molecular  weight. 

The  principal  objective  of  the  present  study  was  to  obtain 
time-average  flow  field  data  in  sufficient  detail  to  define  the 
recirculating  flow  fields  that  develop  from  ducted,  non-reactlve 
but  variable-density,  axisymmetric , turbulent  jet  mixing.  The 
density  variations  were  achieved  by  utilizing  air  as  the  primary 
stream  fluid  and  hydrogen  as  the  secondary  stream  fluid.  The  data 
obtained  were  spatial  distributions  of  time-average  velocity,  hydro- 
gen mass  fraction,  and  wall  static  pressure  for  three  sets  of  pri- 
mary and  secondary  stream  influx  conditions.  The  data  may  be  used 
for  preliminary  design  estimates  for  combustors,  or  combustor  pro- 
cesses, and  to  develop  and  evaluate  theoretical  models  or  predictive 
techniques  for  such  flows  and  processes. 

In  the  present  study,  the  data  obtained  were  compared  with 
theoretical  predictions  made  by  using  a slightly  modified  version 
of  the  theoretical  model  and  numerical  solution  procedure  developed  by 
Gosman,  et  al.,  (Ref.  1).  The  theoretical  solutions  were  calculated 
based  on  a simple  Prandtl-type  eddy  viscosity  model  of  the  turbu- 
lent transport  processes. 

1.2  SURVEY  OF  THE  LITERATURE 

Thrlng  and  Newby  (Ref.  2)  were  concerned  with  analyzing  the 
operating  conditions  within  an  oil-fired,  recirculating-flow  burner. 

The  principal  objective  of  their  analysis  was  "...  to  employ  an 
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Isothermal  constant  density  mixing  model  to  predict  approximately  the 

position  of  a stoichiometric  contour  in  a hot  flame  burning  fluid 

atomized  oil."  The  importance  of  the  paper  by  Thring  and  Newby  was 

that  it  presented  evidence  that,  when  R /R  > 10,  the  jet  mixing 

s p 

processes  in  combustors  could  be  interpreted  qualitatively  like  free 
turbulent  jet  mixing  and  that  they  identified  a parameter  that 
"correlated"  combustor  mixing  processes  with  variations  in  the  ratios 
of  the  weight  flows  and  radii  of  the  primary  stream  to  the  secondary 
stream. 

Curtet  (Ref.  3)  presented  an  analysis  of  ducted,  turbulent  jet 
mixing  processes  which  was  developed  by  himself  and  A.  Craya  based  on 
solving  the  boundary  layer  forms  of  the  Navier-Stokes  equations  using 
an  integral  technique.  In  order  to  obtain  solutions,  Curtet  assumed 
that  the  confined  jet  mixing  processes  could  be  based  on  empiricisms 
developed  by  analysis  of  free  turbulent  jet  mixing  processes.  Curtet 
showed  that  this  type  of  analysis,  an  integral  analysis,  could  predict 
certain  features  of  the  process  of  confined  jet  mixing  with  recircula- 
tion, for  example,  the  amount  of  material  in  the  recirculation  eddy, 
over  appropriate  ranges  of  characteristic  parameters  of  the  system, 
such  as,  Ug  , u , R , and  R . The  analysis  by  Craya  and  Curtet  was 
developed  specifically  for  constant  density,  isothermal,  confined, 
turbulent  jet  mixing.  Curtet,  like  Thring  and  Newby,  reported  a 
mixing  similitude  parameter  that  would  serve  to  correlate  the  gross 
structural  features  of  the  flow  fields  as  they  varied  with  changing 
primary  and  secondary  stream  Influx  conditions. 

To  supplement  the  theoretical  study,  Curtet  presented  the  results 
of  an  experimental  Investigation  of  recirculating  flow  in  a combustor 
configuration.  The  experimental  investigation  covered  recirculating 
flows  in  both  planar  and  axlsymmetric  combustor  configurations.  The 
data  were  obtained  from  experiments  based  on  isothermal,  constant 
density,  confined,  turbulent  jet  mixing  but  unfortunately,  the  data 
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may  be  dependent  on  the  exit  geometry  of  the  combustor  that  was  used 
in  the  experiments  and  thus  may  lack  generality. 

On  the  other  hand,  papers  providing  detailed  data  on  confined  jet 
mixing  with  recirculation  that  do  not  depend  on  downstream  or  exit 
conditions  were  prepared  and  presented  by  Becker,  Hottel,  and  Williams 
(Refs.  4,  5,  and  6).  The  papers  presented  data  which  define  the  wall 
static  pressure,  concentration,  and  velocity  fields  which  result  from 
constant  density,  isothermal,  confined  jet  mixing  processes  with  re- 
circulation. The  data  also  serve  to  define  the  intermittent,  or 
"turbulent,"  nature  of  this  class  of  flows  because  Becker,  Hottel,  and 
Williams  presented  spatial  distribution  of  intermittency  and  concen- 
tration fluctuation  intensities  that  existed  in  their  recirculating 
combustor  flows.  Becker,  Hottel,  and  Williams  proposed  yet  another 
parameter  (C^)  to  serve  as  a parameter  for  the  correlation  of  the 
gross  flow  field  properties  with  variations  in  primary-to-secondary 
stream  velocity  and  radius  ratios.  Their  definition  of  C^,  which  they 
named  the  Craya-Curtet  parameter  (Ref.  4),  is 

U, 


‘k 

—fr 

u 


(2) 


where 


1 2 
-iy[u  R + 
P P 


- Rp] 
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The  form  of  given  above  is  based  on  uniform  parallel  Influx  or 
entrance  conditions  of  the  primary  and  secondary  streams  when  they  are 
separated  by  a thin  nozzle  lip.  Details  and  analysis  of  the  Becker, 
Hottel,  and  Williams  experiments  are  provided  in  Ref.  7. 
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Barchilon  and  Curtet  (Ref.  8)  also  presented  experimental  data 
defining  the  flow  fields  developing  from  isothermal,  constant  density, 
confined,  axisymmetric , turbulent  jet  mixing  with  recirculation.  Their 
data  defined  time-averaged  axial  velocity  fields,  wall  static  pressure 
distributions,  and  turbulent  velocity  fluctuation  Intensity  distributions. 
Barchilon  and  Curtet  also  compared  and  correlated  the  spatial  variations 
of  the  gross  flow  field  characteristics,  such  as  the  structure  of  the 
recirculating  eddy,  with  the  values  of  the  Craya-Curtet  parameter  (C^) . 
Their  stuJy  was  carried  out  with  both  gaseous  and  liquid  streams  as 
the  experimental  fluids. 

Other  papers  have  been  written  about  recirculating  combustor  flows, 

but  these  papers  are  not  sufficiently  useful  to  warrant  discussion  in 

the  present  study.  At  present,  it  is  more  important  to  consider  the 

differences  and  similarities  in  the  data  presented  by  Becker,  Hottel, 

and  Williams  and  by  Barchilon  and  Curtet  since  their  papers  contain  the 

principal  data  which  define  ducted  jet  mixing  with  recirculation  under 

isothermal,  constant  density  conditions.  In  both  cases,  experiments 

were  conducted  in  duct  geometries  where  the  duct-to-pr imary  jet 

radius  was  relatively  large  (R  /R  ^ 10)  so  that,  to  a fair  approxi- 

s p 

mation,  the  jet  mixing  processes  exhibit  qualitative  and  quantitative 
features  similar  to  free  turbulent  jet  mixing  processes.  Both  set  of 
experiments  were  concerned  with  the  investigation  of  constant  density, 
isothermal,  axisymmetric,  ducted  turbulent  jet  mixing  with  recirculation 
wherein  both  primary  and  secondary  streams  enter  the  mixing  region 
with  uniform  velocity  profiles.  Both  sets  of  experiments  provided 
data  defining  the  time-averaged  axial  velocity  fields,  wall  static 
pressure  distributions,  and  distributions  of  quantities  that  characterize 
the  fluctuations  of  properties  of  the  flow  fields  in  an  attempt  to 
characterize  the  "turbulent"  nature  of  the  flows.  In  both  cases,  the 
characteristic  features  of  the  flow  fields,  such  as  the  axial  locations 
of  the  forward  and  rear  stagnation  points,  were  correlated  with  C.^,. 

The  similitude  parameter  was  varied  over  a wide  range,  including  conditions 
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resulting  in  blow-off  of  the  recirculation  zone.  Barchilon  and  Curtet 
showed  that  the  instantaneous  structure  of  recirculating  flows  is  not 
regular,  that  is,  not  characterized  by  a single  eddy.  Rather  the 
recirculation  region  is  characterized  by  one  or  more  short-lived 
eddies  and  extremely  high  fluctuating  velocities  (relative  to  local 
or  mean  velocities),  such  that  large-scale  turbulent  motions  exist 
over  random  time  scales.  Becker,  Hottel,  and  Williams  showed  that 
integral  length  scales  or  standard  deviations  of  turbulent  fluctuations 
of  flow  properties  in  jet  mixing  i.lows  with  recirculation  vary  both 
radially  and  longitudinally  throughout  the  flows.  Moreover,  it  was 
shown  that  these  measures  of  turbulent  fluctuation,  such  as  the 
standard  deviation  of  the  axial  intermitting  profiles,  could  be  the 
same  order  of  magnitude  as  the  duct  radius. 

On  the  other  hand,  an  analysis  and  evaluation  of  the  sets  of  data 
presented  by  Becker,  Hottel,  and  Williams  and  by  Barchilon  and  Curtet 
(Ref.  7)  demonstrated  both  qualitative  and  quantitative  differences 
between  sets  of  data  for  conditions  that  were  supposedly  equivalent, 
based  on  correlations  with  C^.  For  example,  while  the  upstream  edge 
of  the  recirculation  eddy  (Xpg)  appeared  to  correlate  with  C^, 
behaving  approximately  the  same  in  both  sets  of  experiments,  the  data 
defining  the  behavior  or  correlation  of  the  downstream  stagnation 
point  location  (X^g)  with  was  substantially  different  in  the  two 
sets  of  experiments.  Also,  the  rates  of  decay  of  axial  velocity  on 
the  centerline  of  the  ducts  differs  qualitatively  between  the  two  sets 
of  experiments  for  turbulent  mixing  conditions  that  are  supposedly 
equivalent,  based  on  correlations  with  C^.  Therefore,  it  is  clear 
that  additional  experiments  are  required  to  determine  if  mixing 
similitude  parameters  such  as  the  Craya-Curtet  parameter,  are  suitable 
for  expressing  the  equivalency  of  confined,  turbulent  flows.  Also,  the 
data  presented  by  Becker,  Hottel,  and  Williams  and  by  Barchilon  and 
Curtet  which  illustrate  the  "turbulent"  nature  of  recirculating  flow 
fields  must  be  verified  and  extended  by  future  research  on  such  flows. 
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In  practical  applications,  the  flows  in  most  recirculating  flow 
combustors  are  established  by  variable  density  turbulent  jet  mixing. 

The  density  variations  usually  result  from  the  mixing  of  primary  and 
secondary  streams  of  different  temperatures,  molecular  weights,  or 
from  exothermic  chemical  reactions  occurring  in  the  mixing  process. 

There  are,  apparently,  no  experimental  data  defining  variable  density, 
confined,  turbulent  jet  mixing  with  recirculation  which  exist  in 
sufficient  detail  to  be  useful  for  engineering  design  purposes  or  for 
the  development  and  evaluation  of  theoretical  predictive  techniques  for 
combustor  design.  Of  immediate  Interest  is  the  mean  velocity,  con- 
centration, and  wall  static  pressure  field-definition  for  recirculating 
combustor  flows  wherein  large  density  gradients  initially  exist  across 
the  jet  mixing  layers.  Experiments  may  be  either  reactive  or  nonre- 
active, but  the  relative  simplicity  of  nonreactive  to  reactive  experi- 
ments strongly  implied  that  nonreactive  experiments  be  conducted  first. 
Therefore,  a primary  purpose  of  the  present  study  was  to  provide  a 
"first-cut"  set  of  data  defining  the  flow  fields  of  variable  density, 
nonreactive,  ducted,  axisymmetric  turbulent  jet  mixing  with  recirculation 
where  the  density  variations  were  caused  by  the  mixing  of  two  streams 
of  greatly  different  molecular  weight. 

2.0  THEORETICAL  ANALYSIS 


2.1  INTRODUCTION 

The  special  code  for  combustor  problems  developed  and  reported  by 
Gosman,  et  al.  (Ref.  1),  was  designed  to  predict  the  flow  fields  within 
general  axisymmetric  combustor  configurations.  The  code,  called 
hereafter  simply  the  special  code,  was  applied  in  the  present  investi- 
gation together  with  an  assumed  model  for  the  turbulent  or  effective 
viscosity  to  make  theoretical  predictions  of  the  variable-density, 
combustor-type,  recirculating  flow  fields  experimentally  investigated 
in  the  present  study. 
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A set  of  governing  equations  of  motion  for  viscous,  axisymmetric , 
recirculating  flows  was  derived  from  the  Navier-Stokes  equations,  the 
species  conservation  equation,  and  an  identity  relating  stream  function 
and  vorticity  by  Gosman,  et  al.  The  equations  are,  respectively, 


^(r  + rd 


= 0 


(5) 


9 , 1 9 , 1 8,|;. 

^ + ro)  = 0 


(7) 


The  variable  (<(/)  is  the  stream  function,  defined  by  the  set  of  equations 


9tJ; 

■5:3^  = - pvr 


(8) 


4 


and 


9r 


pur 


(9) 


The  variable  (w)  is  the  vorticity  divided  by  the  radius,  that  is. 


0)  * fi/  r 


(10) 
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where  vorticity  was  written,  as 


9v 

9x 


9u 

T? 


(11) 


The  variable  (C)  represents  any  scalar,  "transportable,  conserved" 
quantity.  For  example,  in  the  present  investigation,  it  represents 
the  mass  fraction  of  the  secondary  stream  fluid  which  is  gaseous  hy- 
drogen. 

The  governing  equations  are  a set  of  coupled,  elliptic,  partial 
differential  equations.  Therefore,  solution  can  be  obtained  only  if 
boundary  conditions  for  ip,  w,  and  ? are  prescribed  over  the  entire 
boundary  enveloping  the  flow  field. 


2.2  NUMERICAL  SOLUTION  TECHNIQUE 


A method  of  solving  the  set  of  governing  equations  was  developed 
and  presented  by  Gosman,  et  al.  (Ref.  1).  They  cast  the  set  of  governing 
equations  into  a common  format,  specifically,  writing  each  equation  in 
the  form 


9r^'  97^^ 


lb.r-^(c,<p) 
9x  ‘ (P  9x''  4)^-' 


= 0 


(12) 


where  p represents  0,  u,  or  The  necessary  and  corresponding  set 
of  coefficients  (a^,  b^,  and  c^)  are  those  shown  in  Table  1 opposite 
their  appropriate  functions  (<()).  The  terms  (d^)  are  the  corresponding 
source  terms  for  the  4 and  are  given  in  Table  2 opposite  their 
corresponding  functions  (p). 
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The  source  term  for  o)  contains  a function  (S^)  given  by 

= |•[grad(i^•V)  •isoCi^'grad 

+ grad(i^*V) •isoCi^'grad  y^)-  iso(div  V) ‘grad  y^]  (13) 


where  i^  and  i^  are  unit  vectors  in  the  x and  r direction,  respectively, 
and  where,  given  the  dummy  variables  (A  and  X) 


grad  A = i i 

® 9x  X ar  r 


(14) 


iso  A = -|^i  + |^i 

9r  X 9x  r 


(15) 


and 


3A 


div  A - ifr  ^(rA^)] 


(16) 


in  cylindrical  polar  coordinates  for  axisymmetric  systems. 


Therefore,  in  an  expanded  format  becomes 


S 


u 


■ ax-* 


3u 

rf  - 


3v 


♦ 


3 ,3u 
'5x“^7r^3x 


3 ,3u 


+ 


(17) 
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Gosman,  et  al. , argue  that  contributes  negligibly  to  the  total  vorticity  [ 

source  term  (d  ) and  may  be  neglected.  Certainty,  this  is  true  for  con- 
0) 

stant  viscosity  flows.  In  the  present  study,  was  presumed  to  be  neg- 
ligible despite  the  fact  that  considerable  viscosity  gradients  existed  ; 

in  the  fields  of  the  numerical  solution  because  of  the  turbulence  transport 
model  used  in  the  present  study. 

Gosman,  et  al. , presented  a technique  for  obtaining  a set  of  finite- 
difference  equations  analogous  to  and  derived  from  Eq.  (12).  The  ultimate 
result  of  their  analysis  is  a set  of  coupled,  algebraic  equations  called 
successive  substitution  equations  or  formulas.  One  such  equation  exists 
for  each  dependent  variable  (((>).  The  method  of  solving  for  the  values  of 
(J)  in  the  interior  of  the  computational  region  is  to  apply  the  successive 
substitution  formula  alternately,  in  complete  sweeps  of  the  computational 
field,  in  a point  by  point  manner  for  each  (p.  The  stability  of  the  com- 
putational procedure  is  obtained  for  flows  characterized  by  large  (greater 
than  two)  Reynolds  numbers  that  are  based  on  computational  cell  spacing, 
by  (1)  the  use  of  an  "up-wind"  differencing  technique  for  establishing 
finite  difference  expressions  of  convective  flux  terms  of  Eq.  (12),  (2) 
by  requiring  that  certain  coefficients  in  the  successive  substitution 
formula  be  bounded  in  their  individual  and  collective  numerical  values; 
and,  finally  (3)  by  an  appropriate  use  of  under-  or  over-relaxation 
techniques.  In  the  present  study,  no  systematic  investigation  or  use 
of  under-  or  over-relaxation  factors  was  attempted.  The  under-  or 
over-relaxation  factors  that  were  used  ranged  from  0.9  to  1.1. 

In  addition,  in  the  course  of  the  present  study,  the  numerical 
solutions  of  the  flow  field  were  accepted  as  being  converged  when  the 
maximum  variations  in  the  axial  distribution  of  centerline  axial 
velocity  ceased  to  change  by  more  than  one  percent  with  successive 
iterations  of  the  flow  field. 
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2.3  BOUNDARY  CONDITIONS 

Techniques  for  the  establishment  of  boundary  conditions  for  the 
dependent  variables  (t(/,  w,  and  ?)  are  provided  by  Gosman,  et  al.  In 
general,  at  the  solid  walls,  the  value  of  w and  the  gradients  of  o)  are 
prescribed  such  that  axial  and  radial  velocity  components  are  zero.  At 
the  plane  of  the  primary  jet  nozzle  exit,  the  boundary  conditions  on 
ii  and  0)  are  such  that  they  define  uniform,  parallel  primary  and  secondary 
flow  influx  conditions  everywhere  except  at  the  nozzle  lip.  Along  the 
axis  of  symmetry,  the  special  code  establishes  the  axis  value  of  w as  a 
linear  extrapolation  of  the  values  of  u>  at  one  and  two  radially  adjacent 
nodes.  At  the  outlet  plane,  the  axial  gradients  of  \p  and  u are  set  to 
zero.  For  the  conserved,  transportable  scalar  quantity  (?)  the  boundary 
condition  on  the  wall  and  on  the  axis  of  symmetry  is 

BC/dr  = 0 (18) 


At  the  plane  of  the  primary  nozzle  exit,  C is  constant  across  the  primary 
and  secondary  stream,  equal  to  0.0  and  1.0,  respectively.  Finally  at 
the  outlet  plane,  the  boundary  conditions  are  such  that 

8C/3X  = 0 (19) 


2.4  PRACTICAL  DIFFICULTIES  ARISING  IN  THE  APPLICATION  OF  THE 
SPECIAL  CODE 

Four  major  problems  arose  in  the  application  of  the  special  code 
to  predict  the  turbulent,  recirculating  flows  considered  in  the  present 
study: 

1.  The  specification  of  vorticity  at  a nozzle  lip. 
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2.  The  prediction  of  static  pressure  fields, 

3.  The  numerical  problems  that  arise  when  making 
computations  by  utilizing  highly  stretched  or 
nonuniform  mesh  configurations,  and 

4.  The  problem  of  obtaining  converged  solutions 
for  flows  with  large  density  gradients  in  the 
fields. 

Specific  approaches  have  been  taken  to  resolve  the  difficulties. 

The  approaches  which  are  summarized  below  are  presented  in  detail  in 
Ref.  7. 

2.4.1  Boundary  Condition  for  Vorticity  at  a Protruding  Corner 

Gosman,  et  al.,  gave  techniques  and  methodology  for  establishing  or 
defining  boundary  values  for  the  dependent  variables.  However,  it  was 
found  (Ref.  7)  that  the  vorticity  at  the  lip  of  the  primary  nozzle  was 
one  of  the  principal  parameters  that  determined  the  development  of  the 
theoretical  solution,  that  is,  of  the  development  of  the  downstream  flow 
field.  Gosman,  et  al.,  incorporated  a special  procedure  to  specify  w at 
the  edges  of  the  nozzle  because,  as  they  noted,  this  point  is  a point 
of  flow  separation  from  a "protruding  corner." 

Figure  2 shows  the  geometry  of  the  boundary  of  the  flows  under 
investigation,  and  as  noted,  the  points  CA  and  CB,  at  the  edges  of 
the  nozzle  lip,  define  the  points  of  flow  separation.  Many  techniques 
exist  for  the  specification  of  vorticity  at  such  corners.  For  example, 
Roache  (Ref.  9)  lists  seven  schemes  for  assigning  values  to  o)  at  such 
corners.  Gosman,  et  al. , sets  the  value  of  the  stream  function  at  the 
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points  CA  + 1 and  CB  + 1 equal'  to  the  values  of  stream  function  at  CA 
and  CB,  respectively; 


CA 


(20) 


and 


’^CB  + 1 “ '^CB 


Then,  the  value  of  the  corner  vorticity  at  CA  and  CB  is  established 
through  a boundary  function  that  lepresents  the  condition 


“CA  " ° ^CA  ~ '^CB  " 


(21) 


(22) 


However,  the  points  (CA  + 1 and  CB  +1)  now  become  effective  corner 
points.  The  value  of  the  vorticity  at  these  points  is  established 
through  the  following  special  routine.  In  principle,  Eq.  (7)  applies, 
but  in  reverse. 


r.o.  = - [ 9 (J-  . 9 (-L  |i)  1 

'■  3x'-pr  Jx-'  9r'-pr  Tr-*  ^ 


(23) 


Therefore,  Gosman  et  al.,  evaluate  this  equation  numerically  and 
approximately,  for  example,  at  CA  + 1 , by 


(jj 


ip  \p  ^ip  ^\p 


(24) 


where  Cg  refers  to  the  coefficient  of  the  contribution  to  the  stream 

'I' 

function  equation  (Eq.  (7))  from  the  eastern  node  (which  in  the  present 
special  routine  is  CA) . The  other  coefficients  have  corresponding 
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meanings.  Thus,  the  "corner"  vorticity  has  been  transferred  to  the 
first  node  point  downstream  at  the  same  radial  position.  As  a result, 
the  jet  must  issue  from  the  nozzle  without  expanding  or  contracting  up 
to  the  first  node  point  directly  downstream  of  the  nozzle  lip.  But,  the 
mathematical  relations  thus  established  by  Eqs.  (20),  (21),  and  (23)  make 
corner  or  nozzle  lip  vorticity  dependent  on  the  computational  mesh 
configuration  used  to  simulate  the  nozzle  geometry,  the  thickness  of 
the  simulated  nozzle  lip,  and  the  x-wise  mesh  spacing  separating  such 
points  as  CA  and  CA  + 1 . Such  difficulties  have  also  been  noted  by 
other  investigators,  for  example,  by  Saccia  and  Kennedy  (Ref.  10), 
without  being  resolved.  Investigation  (Ref.  7)  has  shown,  however, 
that  the  establishment  of  a representative  corner  or  lip  vorticity  is 
very  important  and,  at  present,  must  be  established  by  trial  and  error. 


Therefore,  in  the  present  study,  an  approach  was  taken  which 
attempted  to  relate  corner  vorticity  to  the  slope  of  the  primary 
nozzle  velocity  profile,  that  is,  at  the  nozzle  lip. 


’"“CA  " ^CA 


= X u /R 
P P 


9u 

9r 


nozzle 


lip 


(25) 


Hopefully,  therefore,  the  constant  (X)  can  be  obtained  by  inspection 
of  the  primary  flow  nozzle  exit  velocity  profiles,  which  are  experimentally 
determined . 

2.4.2  Computation  of  Static  Pressure  Fields 

The  method  used  in  the  present  study  to  compute  static  pressure 
fields  has  been  applied  to  both  constant  density  and  variable  density 
ducted  jet  mixing  flows  with  recirculation  (Ref.  7).  This  method  is 
different  from  the  method  Incorporated  in  the  special  code  wherein  the 
static  pressure  field  is  obtained  by  integrating  the  Navler-Stokes 
equations  once  the  velocity,  density,  and  effective  viscosity  fields 
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have  been  established.  Instead,  the  method  used  in  the  present  study 
for  constant-area  ducted  flows  was  found  to  be  much  simpler  in  procedure 
and  at  least  as  accurate  as  the  method  of  Gosman,  et  al.  The  method  is 
based  on  two  assumptions: 

1.  Static  pressure  is  constant  across  the  duct,  and, 
hence  is  a function  only  of  axial  distance,  and 

2.  Wall  shear  stresses  contribute  negligibly  to  the 
sum  of  pressure  area  forces  and  momentum  flux  in 
an  application  of  Newton's  Second  Law  to  the  duct 
flow. 

Under  these  conditions,  the  impulse  function  can  be  evaluated  at  any 
station  of  the  duct,  and  set  equal  to  its  value  at  the  inlet  station, 
that  is, 

F(x)  = F(o)  (26) 

where  F is  the  impulse  function  defined  as 

F(x)  = / (p  + pu^)  dA  (27) 

A(x) 


Applying  Eq.  (6)  at  the  inlet  station  (x  = 0)  and  any  other  axial 
station  (x  = L)  and  equating  the  expressions  yield 


u^A 

^ A^/  ^ “^A?rf 

A(L) 


p U ^A 
s s s 

o (28) 

A(L) 
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or 

2„2  p R 

puR  ss  , , 7 r ^ 7 

p(L}  = p(o]  + P P P + p (Rf  - - — I pu  rdr  (29) 

2 2 P 2 1 

R'^  R*^  R'^  J 

s s s "o 

Thus,  the  pressure  at  any  station  of  the  duct  can  be  found  by  integration 

2 

of  the  profile  of  pu  at  that  station  in  the  duct. 

2.4.3  Computations  Using  Nonuniform  Mesh  Networks 

The  flow  fields  to  be  predicted  are  formed  by  turbulent  jet  mixing 
between  a relatively  small  diameter,  axisymmetr ic  jet  and  a relatively 
much  larger  diameter  ducted  secondary  stream.  It  was  necessary,  in 
order  to  have  reasonable  computation  times,  to  carry  out  the  numerical 
solution  in  computational  mesh  networks  that  were  of  nonuniform  node 
spacing  (Ref.  7).  The  region  near  the  primary  jet  nozzle  exit  plane, 
called  the  near  field  of  the  jet  mixing  region,  was  defined  by  a mesh 
network  with  node  spacings  which  were  relatively  small  compared  to  the 
node  spacings  out  near  the  duct  wall  and  downstream  in  zones  where 
recirculation  was  expected  to  occur.  Unfortunately,  computations  made 
utilizing  highly  nonuniform  networks  of  arbitrary  construction  may  result 
in  poor  quality  solutions  (Ref.  7).  For  example,  in  computations  of 
nonreactive,  multispecies  flows,  it  has  been  observed  (Ref.  7)  that  the 
conservation  of  species  criterion  may  be  violated  by  the  solution  pro- 
cedure; the  species  concentrations  at  the  nodes  converge  to  values  which 
fail  to  satisfy  requirements  of  global  conservation.  However,  such 
errors  can  be  lessened  by  trial-and-error  manipulation  of  the  mesh 
network  configuration  in  which  the  solution  is  obtained.  The  errors 
can  be  lessened  to  the  extent  that  approximate  species  conservation  can 
be  obtained  using  the  solution  procedure.  For  example,  mesh  configura- 
tions were  developed  (Ref.  7)  for  the  computation  of  the  flows  reported 
by  Becker,  Hottel,  and  Williams  and  for  the  variable-density  flows 
reported  herein. 
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2.4.4  Obtaining  Converged  Variable  Density  Solutions 

When  attempts  were  made  to  compute  the  hydrogen-air  recirculating 
flows,  it  was  found  that  the  numerical  solution  procedure  would  not 
converge.  The  difficulty  was  attributed  to  numerical  Instabilities 
that  arose  because  of  the  arbitrary  initial  distributions  of  the  stream 
function  and  mixture  fraction  (or  density)  that  were  used  to  start  the 
calculation  procedure.  In  order  to  circumvent  this  problem,  the  "density 
dependence"  of  the  stream  function  was  eliminated  in  the  beginning  of 
the  iteration  process  (Ref.  7).  To  implement  this  procedure,  the  stream 
function  variable  was  modified  so  that  it  was  explicity  independent  of 
density  in  the  first  several  hundred  iterations  of  the  solution  procedure. 
The  net  result  of  the  modification  was  to  impose,  initially,  an  incorrect 
global  continuity  equation  which  linked  the  density  field  to  the  field 
of  the  modified  stream  function  in  a very  indirect  or  weak  relationship. 
This  resulted  in  initially  stable  numerical  solution  procedures,  albeit 
with  improper  governing  equations.  After  400,  or  so.  Iterations  an 
approximately  realistic  density  and  velocity  field  was  obtained,  and 
at  this  point  in  the  solution  procedure,  integration  of  the  pu  field 
yielded  an  initial  stream  function  field.  This  field  of  density-dependent 
stream  function  was  sufficiently  well  posed  to  allow  continued  convergence 
of  the  solution  procedure  with  the  correct  governing  equations.  Details 
of  this  method  for  securing  convergence  of  the  solution  procedure  in 
variable  density  computations  are  outlined  in  Ref.  7. 

2.5  TURBULENCE  TRANSPORT  MODEL 

The  closure  problem  which  underlies  the  overall  problem  of  pre- 
dicting turbulent  flows  is  that  the  turbulent  transport  of  momentum, 
energy,  and  chemical  species  must  be  related  to  time-averaged  flow  field 
quantities  or  to  their  gradients.  There  are  several  valuable  references 
which  present  work  done  to  model  turbulence  transport  phenomena  for  a 
wide  variety  of  flows,  (Refs.  11  through  18).  However,  comparatively 
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little  work  has  been  performed  on  turbulence  transport  modeling  in 
recirculating  combustor  flows  other  than  the  simple  observations  that 
the  mixing  is  somewhat  like  the  mixing  of  free  turbulent  jets.  Therefore, 
the  approach  taken  in  the  present  study  was  similar  to  that  proposed  in 
Refs.  3 and  19  and  was  developed  in  Ref.  7.  Two  models  for  effective 
viscosity  were  defined  that  were  based  on  a Prandtl-type  eddy  viscosity 
model  for  free  turbulent  jet  mixing  processes.  One  model,  called  the 
model  without  potential  core,  is  defined  such  that 

Uc(x,r)  = 0.0285  p(x,r)  u R 

n P P 


core 


, X < X 


core 


(30) 


= 0.0285  p(x,r)UpRp(l  + 0 . 0036 /Rp) , x > X^^^^ 

(31) 

Here,  represents  a length  characteristic  of  a potential  core  of  a 

free  turbulent  jet,  namely,  X = 12  R . Note  that  p_  is  constant 
across  the  flow  at  a given  axial  station  in  constant  density  flow.  The 
second  model,  called  the  model  with  potential  core,  is  defined  such  that 

Pc(x,r)  = 0.0285  p(x,r)  u R^  • — , x < X 

p p X — core 

core 

yg(x,r)  = 0.0285  p (x , r) UpRp ( 1+0 . 0036 (X- /Rp ) , x > X^^^^ 

(33) 

in  the  "turbulent"  regions  of  the  flow  field,  and 

Ug(x,r)  = 

in  the  "laminar"  regions  of  the  flowfleld.  The  "laminar"  and  "turbulent" 
regions  of  the  flow  field  are  defined  in  Fig.  3.  The  purpose  of  the 
model  with  potential  core  was  to  investigate  the  Importance  of  accounting 
for  the  non-turbulent  regions  of  the  recirculating  flow  (Ref.  7). 
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Unfortunately,  the  model  with  potential  core  requires  more  empirical 
information  in  its  formulation  than  the  model  without  potential  core 
since  the  spatial  variation  of  the  laminar  regions  of  recirculating 
flows  must  be  specified. 


3.0  APPARATUS 


3.1  TEST  CELL 


To  obtain  data  on  variable  density  mixing  with  recirculation, 
specifically  the  hydrogen-air  mixing  flows,  a combustor  simulator  was 
designed  and  constructed.  The  combustor  design,  shown  in  Fig.  4, 
consists  of  a 5-ft  steel  duct  with  an  inside  diameter  of  5.24  in.  A 
mechanically  driven,  axially  traversing  nozzle  assembly,  shaped  like 
a cylindrical  piston,  was  mounceJ  inside  the  duct.  The  nozzle  assembly 
(Fig.  4)  was  traversed  in  discrete  steps  during  the  testing  to  vary  the 
distance  (L)  between  the  primary  jet  lozzle  exit  plane  and  the  radially 
positionable  pitot  pressure  and  gas  sampling  probe.  An  0-ring  seal 
system  isolates  the  flow  in  the. duct  from  the  atmosphere. 


The  nozzle  assembly  includes  a circular,  0. 524-ln. -diam  central 

nozzle,  which  introduces  the  high  subsonic  velocity,  primary  air  jet, 

and  an  annular  secondary  injector  assembly  through  which  the  outer 

hydrogen  stream  enters  the  duct.  The  duct-to-jet  radius  ratio  (R  /R  ) 

s p 

is  10,  Details  of  the  primary  flow  nozzle  are  shown  in  Fig,  5.  In 

Fig.  4,  it  is  seen  that  the  secondary  flow  must  pass  through  two  porous 

plates  and  a screen.  The  two  porous  plates  were  Installed  to  reduce 
the  total  pressure  of  the  secondary  stream  to  appropriate  levels.  The 
screen  with  57-percent  porosity  was  installed  to  distribute  the 
secondary  mass  flow  uniformly  at  the  primary  nozzle  exit  station. 
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3.2  INSTRUMENTATION 

3.2.1  Pressure  and  Temperature  Measurement 

A pitot  pressure  and  gas  sampling  probe  having  both  upstream  and 
downstream  facing  probe  tips  was  installed  in  the  duct  at  a fixed 
axial  position.  The  probe  can  be  traversed  radially  across  the  duct 
by  a geared  drive  mechanism.  The  probe  dimensions  are  shown  in 
Fig.  6.  By  combining  the  motion  of  the  axial-traversing  nozzle 
assembly  and  the  radially  traversing  probe,  axial  and  radial  surveys 
were  made  of  the  jet  mixing  and  recirculating  flow  fields.  In 
addition  to  the  pitot  probe  assembly,  a total  of  50  static  pressure 
orifices  were  installed  along  the  duct  wall,  the  majority  of  which 
were  spaced  about  one  inch  apart,  in  a helical  pattern  that  completes 
one  revolution  of  the  duct  in  a distance  of  45.5  in.  Orifices  are 
spaced  closer  than  one  inch  in  the  vicinity  of  the  probes.  Both  static 
pressure  and  gas  samples  were  obtained  through  the  wall  orifices.  The 
static  pressure  orifices  were  connected  to  a water-filled  manometer 
board  which  was  photographed  during  the  tests.  When  gas  composition 
measurements  were  required  at  the  wall,  a scanner  valve  system  connected 
the  wall  orifices  to  the  gas  concentration  measurement  system.  Additional 
wall  orifices  were  in  the  planes  of  the  tips  of  the  pitot  probes,  and 
measurements  from  these  were  utilized  in  the  calculation  of  average  axial 
velocity . 

Probe  pitot  (total)  pressures  were  measured  using  a system  composed 
of  two  transducers.  In  the  pressure  differential  range  from  0 to  10  mm 
of  mercury,  the  pitot  pressures  were  sensed  with  a variable  capacitance 
transducer,  while  in  the  range  from  10  to  215  mm  of  mercury  the  pitot 
pressures  were  sensed  by  a 0-  to  4-psi  strain-gage  transducer.  Because 
of  the  sensitivity  of  the  transducer  used  to  measure  the  lower  range  of 
the  pitot  pressure,  axial  velocities  lower  than  100  ft/sec  in  pure  hydrogen 
or  velocities  less  than  25  ft/sec  in  pure  air  could  not  be  accurately 
measured. 
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Temperatures  were  measured  using  copper-constantan  thermocouples. 

The  temperatures  recorded  during  testing  were  the  total  temperature  of 
the  gaseous  supplies  of  primary  and  secondary  fluid,  and  they  were 
measured  upstream  of  critical-flow  venturis,  which  were  used  to  meter 
the  primary  and  secondary  mass  flows.  The  venturis  were  calibrated  so 
that  their  discharge  coefficients  were  known  over  the  ranges  of  metered 
mass  flow  encountered  during  the  test  program. 

3.2.2  Gas  Concentration  Measurements 

Concentrations  of  hydrogen  in  gas  samples  were  measured  using  a 
system  of  gas-sample-driven,  fluidic  oscillators.  The  theory  and  use 
of  fluidic  oscillator  systems  to  measure  gas  concentrations  is  described 
in  Refs.  20,  2l,  and  22.  The  output  frequency  of  a fluidic  oscillator 
depends  on  the  temperature  and  molecular  weight  of  the  gas  sample  passing 
through  the  oscillator,  the  specific  construction  of  the  oscillator,  and 
che  pressure  drop  across  the  oscillator.  However,  through  calibration, 
the  oscillator  output  frequency  can  be  established  as  only  a function  of 
the  molecular  weight  of  the  gas  sample  passing  through  the  oscillator. 
Calibration  of  the  fluidic  oscillator  system  was  accomplished  using  gas 
samples  of  air,  hydrogen,  and  mixtures  of  nitrogen  and  hydrogen  having 
molecular  weights  of  20  and  10. 

3.3  TEST  PROCEDURE 

The  pressure  and  concentration  instrumentation  was  calibrated  before 
each  test  by  applying  to  each  system  a set  of  known  pressures  and  con- 
centrations and  recording  the  corresponding  instrumentation  output  values 
from  each  system.  The  sensitivity  of  each  measurement  system  was  obtained 
over  ranges  of  values  expected  during  testing. 

The  calibration  data  and  the  data  taken  during  each  test  were 
recorded  on  seven-track.  Incremental,  magnetic  tape  by  a programmable 
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Hewlett-Packard  2010  data  acquisition  system  (Ref.  23).  In  special  (■ 

Instances,  in  order  to  obtain  strictly  qualitative  information,  probe  | 

pitot  pressures  and  selected  wall  static  pressures  (sensed  by  0 to  10  mm  | 

of  mer_ury  variable  capacitance  transducers)  were  displayed  on  an  X-Y  i 

recorder. 

I 

‘ 

At  the  beginning  of  a test,  the  desired  primary  and  secondary  mass  I 

flow  rates  were  established  by  setting  and  monitoring  appropriate  air  ] 

' and  hydrogen  gas  supply  pressures  and  temperatures  upstream  of  the  | 

critical-flow  venturis.  The  piston-like  nozzle  assembly  was  positioned  | 

at  a desired  distance  from  the  probe.  Data  were  recorded  as  the  probe  j 

assembly  was  traversed  across  the  flow  in  a series  of  discrete  steps. 

The  nozzle  assembly  was  then  moved  1.38  in.  toward  the  probe,  which  ■ 

corresponds  to  the  distance  separating  the  forward  and  reverse  facing 

pitot  probes  (Fig.  6),  and  another  radial  survey  was  conducted.  Thus, 

forward  and  aft  facing  probes  surveyed  a radial  path  at  the  same  axial 

location  relative  to  the  nozzle  assembly.  Gas  concentrations  were 

sampled  with  the  forward  facing  probe,  and  pitot  pressures  were  taken  , 

with  both  forward  and  aft  facing  probes.  The  complete  pitot  pressure 
and  gas  concentration  profiles  at  a given  axial  station,  relative  to  the 
nozzle,  were  constructed  from  the  surveys  of  the  forward  and  aft  facing 
1 probe  tips  at  this  station. 

At  the  Initiation  of  the  hydrogen-air  mixing  experiments,  radial 
ij  surveys  of  pitot  pressure  and  gas  concentration  were  made  across  the 

duct  from  wall  to  wall  to  ensure  that  buoyancy  forces  would  have  a 
negligible  effect  on  the  mixing  processes.  The  initial  surveys  con- 
firmed that  the  flow  fields  were  symmetrical  about  the  centerline, 
and,  for  the  tests  reported  herein,  radial  surveys  of  int  flow  fields 
were  limited  to  -0.15  ^ 1.  0'®^-  This  procedure  minimize  hydrogen 

consumption  and  test  time. 


f 
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3.4  DETERMINATION  OF  THE  PRIMARY  JET  NOZZLE  FLOW  CHARACTERISTICS 

Before  experimentation  was  initiated,  a pitot  pressure  survey  of 
the  jet  issuing  from  the  primary  nozzle  was  made  to  establish  the  velocity 
field  at  the  inlet  of  the  combustor,  hence,  to  ensure  that  the  primary 
jet  velocity  profile  was  symmetrical.  Typical  velocity  profiles  obtained 
in  orthogonal  surveys  are  shown  in  Fig.  7 and  compared  with  a fully 
developed  turbulent  pipe  flow  profile  at  an  equivalent  Reynolds  number. 

3.5  PRECISION  OF  MEASUREMENTS 

The  uncertainty  of  a measured  parameter  was  determined  from  the 
bias  limits  and  precision  indices  of  the  instrumentation  by  the  proce- 
dures developed  in  Ref.  24.  Given  the  bias  limits  and  precision  indices 
associated  with  the  measurement  necessary  to  define  any  parameter  (<())  the 
total  uncertainty  in  ip  is  defined  as 

^<P  " ^ ^95  (35) 
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is  the  precision  index  for  (ji,  and  where  are  the  measured  parameters 
which  define  ((i;  that  is, 

<P  = (fi(Ci)  , i = 1,  . . . ,N  (38) 

N is  the  number  of  measured  parameters  that  define  (p.  The  parameter 
tg^  is  the  "95  percentile  point  for  the  two-tailed  students  "t" 
distribution." 

The  bias  limits  and  precision  indices  of  Instruments  vary  with 
variations  in  the  values  of  measured  quantities.  Therefore,  represen- 
tative test  conditions  were  selected  upon  which  to  base  the  precision 
indices  and  bias  limits.  Typical  bias  limits,  precision  indices,  and 
uncertainties  in  the  fundamental,  measured  parameters  and  the  axial 
velocity  and  hydrogen  gas  concentrations,  which  are  calculated  from  the 
fundamental,  measured  parameters,  are  presented  in  Table  3. 

4.0  RESULTS  AND  DISCUSSION 

4.1  INTRODUCTION 

The  range  of  test  conditions  established  in  the  hydrogen-air 
experiments  is  presented  in  Table  A.  The  conditions  are  taken  as 
defining  three  sets  of  flows  of  approximately  the  same  velocity  ratio 
, (u  = Ug  three  sets  of  experiments,  the  data  defining 

the  axial  velocity,  hydrogen  concentration,  and  wall  static  pressure 
fields  are  presented  in  Appendix  A.  The  data  have  been  shown  (Ref.  7) 
to  possess  both  differences  from,  and  similarities  to,  the  constant 
density,  recirculating  flow  mixing  data  of  Becker,  Hottel,  and  Williams, 
and  Barchllon  and  Curtet.  In  the  present  study,  experimental  profiles 
of  axial  velocity,  hydrogen  mass  fraction,  wall  static  pressure,  and 
the  stagnation  points  on  the  wall  are  compared  with  corresponding 
theoretical  calculations.  The  theoretical  calculations  were  made  by 
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using  the  theoretical  approach  summarized  in  Section  2.0,  after  applying 
both  viscosity  models;  that  is,  calculations  were  made  with  and  without 
potential  core  simulation. 

4.2  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  PROFILES 
OF  AXIAL  VELOCITY 

The  comparisons  between  experimental  and  theoretical  radial  profiles 
of  axial  velocity  are  shown  in  Fig.  8 for  the  three  sets  of  experiments. 

In  general,  the  shapes  of  the  calculated  velocity  profiles  are  in  agree- 
ment with  experiment  and  show  a jet-like  nature,  as  would  be  expected. 

The  largest  differences  between  calculated  and  experimental  velocity 
profiles  occur  near  the  centerline  of  the  duct.  The  comparisons  between 
calculated  and  experimental  velocity  profiles  do  not  demonstrate  superi- 
ority of  one  viscosity  model  over  the  other.  At  some  axial  stations, 
the  predictions  made  by  using  the  viscosity  model  without  potential  core 
appear  to  be  in  closer  agreement  with  experiment  than  the  calculations 
made  with  potential  core  simulation.  Conversely,  at  other  axial  stations, 
the  reverse  trends  are  observed.  This  behavior  is  noted  in  all  three 
sets  of  comparisons. 

To  obtain  the  calculations  made  using  the  viscosity  model  with  the 
potential  core  simulation,  a considerable  sacrifice  was  incurred  in 
the  rigor  of  the  solution  procedure.  To  obtain  convergent  solutions,  the 
complete  source  term  in  the  vorticity  transport  equation  (Eq.  (12))  had 
to  be  set  to  zero,  i.e., 

d = 0 (39) 

0) 


Moreover,  the  value  of  X in  the  formulation  of  the  corner  vorticity 
(Eq.  (25))  was  reduced  by  approximately  3 percent  from  the  value  used  to 
make  the  calculations  with  a viscosity  model  without  potential  core 
simulation.  Thus,  X changed  from  X = 3.25  to  X = 3.15  when  the  viscosity 
models  were  changed  from  no  potential  core  simulation  to  potential  core 


I 
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simulation.  It  was  also  determined  that  the  values  of  X required  to  I 

make  the  theoretical  calculations  could  not  be  determined  by  inspection 
of  the  primary  jet  exhaust  velocity  profile  as  proposed  in  Section  2.0. 

Instead,  values  of  X were  obtained  by  trial-and-error  calculation  and 
comparison  with  experiment.  The  theoretical  velocity  profiles  shown 
in  Fig.  8 are,  therefore,  correlations  of  experiment  rather  than  pre- 
dictions since  the  values  of  X necessary  to  make  the  computations  could 
not  be  established  beforehand. 

4.3  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RADIAL  PROFILES 

OF  HYDROGEN  MASS  FRACTION  | 

Radial  distributions  of  hydrogen  mass  fraction  obtained  experimentally 
are  compared  with  theoretical  calculations  in  Fig.  9.  The  corresponding 
experimental  hydrogen  mass  fractions  obtained  at  the  wall  have  been 
added  to  this  figure.  The  experimental  and  theoretical  distributions 
are  not  in  satisfactory  agreement,  quantitatively,  but  do  exhibit  similar 
trends.  For  example,  both  the  experimental  and  theoretical  profiles 
indicate  that  the  hydrogen  mass  fraction  tends  to  reach  a maximum  in  the 
Interior  of  the  field,  not  at  the  wall. 

Where  a comparison  of  data  from  tests  C and  D exists  (Fig.  9b),  a 
consistent  difference  is  indicated  by  experimental  data.  The  difference 
in  the  experimental  data  was  determined  to  be  caused  by  a calibration 
error  in  pre-  and  posttest  calibrations  of  the  pitot-gas  sampling 
system  for  test  D.  The  calibration  error  was  accounted  for  in  a 
consistent  manner  in  reducing  data  from  test  D,  but  despite  the  correc- 
tion procedure,  the  difference  between  tests  C and  D persists. 

The  radial  profiles  of  hydrogen  mass  fraction  decay  smoothly 
into  flat,  uniform  profiles  for  all  velocity  ratios.  The  theoretical 
concentration  profiles  show  greater  rates  of  decay,  diffusion,  and 
spreading  than  experiment.  The  calculations  made  using  the  viscosity 
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model  with  potential  core  simulation  appear  to  agree  somewhat  closer 
with  experiment  in  the  near  field,  but  are  in  worse  agreement  with 
experiment  in  the  far  field,  when  compared  with  the  calculations  made 
without  potential  core  simulation.  For  both  viscosity  models,  the 
theoretical  calculations  shown  in  Fig.  9 were  made  with  a numerical, 
effective  Schmidt  number  (Table  1)  of 

(40) 

However,  numerical  experimentation  (Ref.  7)  showed  that  the  theoretical 

predictions  are  not  se.isitive  to  variations  of  up  to  42  percent  in  the 

value  of  the  turbulent  Schmidt  number,  N, 

^c~ 

4.4  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  AXIAL 
DISTRIBUTIONS  OF  AXIAL  VELOCITY  ON  THE  CENTERLINE 


Axial  distributions  of  time-averaged  axial  velocity  on  the  centerline 
of  the  combustor  simulator  are  presented  in  Fig.  10  for  each  of  the 
characteristic  velocity  ratios  (u).  The  experimental  data  exhibit  an 
unusual  feature,  namely,  the  existence  of  regions  in  the  velocity  field 
where  relatively  little  change  occurs  in  the  centerline  axial  velocity, 
for  u = 0.036  and  0.07.  This  feature  appears  as  "bumps"  in  experimental 
centerline  axial  velocity  decay  curves.  The  constant  density  data 
reported  by  Becker,  Hottel,  and  Williams  (Ref.  4)  exhibit  centerline 
axial  velocity  decay  characteristics  somewhat  similar  to  those  in 
Fig.  10  obtained  in  the  variable-density  mixing  experiments.  On  the 
other  hand,  the  theoretical  predictions  show  smooth  decay  for  both 
viscosity  models  for  all  velocity  ratios.  The  reasons  for  this  nature 
of  the  experimental  velocity  fields  are  not  known  although  it  was 
shown  in  Ref.  7 that  the  distributions  of  dynamic  pressure  along  the 
axis  of  the  duct  also  decay  smoothly  for  all  velocity  ratios.  One 
satisfactory  feature  of  the  theoretical  solutions  is  that  they  appear  to 
decay  to  values  in  Cx.yse  agreement  with  experiment  at  large  axial  dis- 
tances from  the  primary  nozzle  exit  plane,  for  example,  at  ^ 120. 
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4.5  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  AXIAL  DISTRIBUTIONS 
OF  HYDROGEN  MASS  FRACTION  ON  THE  WALL  AND  ON  THE  CENTERLINE 

Comparisons  of  theoretical  and  experimental  axial  distributions  of 
hydrogen  mass  fraction  on  the  wall  and  on  the  centerline  are  presented 
in  Fig.  11.  The  experimental  data  from  tests  E and  F show  considerable 
scatter  at  the  wall  near  the  primary  jet  exit  plane.  The  reason  for 
this  scatter  is  not  known  but  experimental  investigation  (Ref.  7) 
eliminated  the  possibility  of  air  leaks  around  the  nozzle  assembly.  An 
interesting  feature  of  the  experimental  data  is  that  the  concentration 
potential  core  extends  further  downstream  than  the  velocity  potential 
core.  The  concentration  potential  core  is  defined  as  the  length  from 
the  exit  plane  of  the  primary  Jet  nozzle  to  the  point  on  the  axis 
of  the  system  where  hydrogen  concentrations  are  first  sensed.  The 
concentration  potential  core  appeared  to  be  about  20  primary  jet 
radii  in  length  (Fig.  11),  whereas  the  velocity  potential  core  appeared 
to  be  about  10  to  12  primary  jet  radii  in  length  (Fig.  10). 

The  theoretical  calculations  of  the  wall  and  centerline  hydrogen 
mass  fraction  are  in  poor  agreement  with  the  experimental  data  for  both 
viscosity  models  in  the  near  field.  Again,  as  in  the  comparisons 
between  experimental  and  theoretical  interior  flow  field  hydrogen  dis- 
tributions, the  theoretical  computations  of  wall  concentration  show 
much  greater  spreading  and  diffusion  of  hydrogen  than  is  observed  ex- 
perimentally. Note,  however,  that  the  theoretical  calculations  of  wall 
concentration  made  using  the  viscosity  model  without  potential  core 
show  good  agreement  with  the  experimental  data  in  the  far  field,  where 
the  radial  concentration  profiles  have  flattened  out.  However,  the 
calculated  wall  concentrations  from  the  viscosity  model  that  simulated 
the  potential  core  regions  clearly  shows  evidence  of  nonconservation 
of  hydrogen  mass  fraction.  Therefore,  further  use  of  this  viscosity 
model  will  require  additional  development  of  computational  mesh  net- 
works as  outlined  in  Section  2.4.3. 
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4.6  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  AXIAL  DISTRIBUTIONS 
OF  STATIC  PRESSURE  ALONG  THE  WALL 


The  experimenta] ly  determined  static  pressure  distributions  are 
compared  with  theoretical  calculations,  made  using  both  viscosity 
models,  in  Fig.  12.  Also  shown  in  this  figure  are  the  maximum  static 
pressures,  or  recovery  pressures,  that  correspond  to  one-dimensional 
flow  theory,  neglecting,  however,  wall  shear  stresses.  The  maximum 
wall  static  pressures  were  calculated  based  on  the  experimental  pres- 
sure and  momentum  flux  conditions  for  each  characteristic  u. 


The  pressure  distributions  computed  with  the  two  viscosity  models 
do  not  compare  well  with  one  another  nor  with  the  experimental  data  in 
the  near  field.  The  calculated  minimum  pressures  in  near  field  pressure 
distributions  increase  as  u decreases  in  the  computations  made  using  a 
viscosity  model  without  potential  core  modeling.  But,  the  calculated 
minimum  pressures  in  the  near  field  static  pressure  distributions  increase 
with  increasing  u in  the  computations  made  with  the  viscosity  model  that 
attempts  to  account  for  laminar  flow  regions.  The  reason  for  this  be- 
havior is  unknown.  In  the  far  field,  however,  where  the  static  pressure 
approaches  its  maximum  value,  both  sets  of  computations  are  in  reasonable 
agreement  with  experimental  data. 


4.7  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  AXIAL  LOCATIONS 
OF  STAGNATION  POINTS  ON  THE  WALL 


Experimental  data  defining  the  axial  velocity  fields  were  inter- 
polated to  determine  points  of  zero  velocity  in  the  three  sets  of  flow 
fields.  The  locations  of  the  zero  velocity  points  in  the  flow  fields 
are  shown  in  Fig.  13.  Because  the  points  of  zero  velocity  are  difficult 
to  determine  experimentally,  these  data  are  the  least  accurate  of  all 
the  experimental  data  presented  herein.  However,  by  using  the  data  in 
Fig.  13,  stagnation  points  on  the  wall  were  defined  and  are  shown  as 
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functions  of  u in  Fig.  14.  Also  shown  in  Fig.  14  are  theoretical  pre- 
dictions of  the  stagnation  point  locations  made  using  both  viscosity 
models.  The  predictions  of  the  forward  stagnation  points  are  in  poor 
agreement  with  the  data,  quantitatively,  but  show  the  same  qualitative 
trends,  namely,  the  stagnation  point  location  moves  downstream  as  u 
Increases.  The  predicted  rear  stagnation  point  locations  are  in  fair 
agreement  with  the  data;  however,  the  inaccuracy  of  the  experimental 
data  defining  the  rear  stagnation  point  locations  is  too  large  to 
confirm  whether  or  not  rear  stagnation  point  varies  with  u as  pre- 
dicted by  the  theory. 


Investigations  (Refs.  4,  5,  and  8)  of  constant  density  mixing  flows 
have  shown  that  the  location  of  the  forward  stagnation  point  varies 
with  the  parameter  C.j,  as  shown  in  Fig.  15.  On  the  other  hand,  the 
data  (Refs.  4,  5,  and  8)  have  not  demonstrated  the  variation  of  the 
rear  stagnation  point  location  with  conclusively  (see  Fig.  15).  To 
be  consistent  with  the  published  results  of  the  constant-density  experi- 
ments, a variable-density  Craya-Curtet  parameter  (C*)  was  defined, 
consistent  with  the  definition  provided  by  Becker,  Hottel,  and  Williams, 
as 
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In  this  formulation,  p*  and  p*  represent  mass  and  momentum  averaged 


densities  based  on  inlet  conditions.  In  particular. 
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However,  it  should  be  pointed  out  that  alternate  definitions  for  p* 

^ and  p*  have  been  proposed  (Ref.  25).  Nevertheless,  by  using  the  above 

definitions  of  p*  and  p*,  values  of  C*  were  determined  and  are  listed 

for  all  of  the  present  experimental  conditions  in  Table  4.  The  variation 

of  the  forward  and  rear  stagnation  point  locations  with  C*  are  compared 

with  the  results  from  the  constant  density  experiments  in  Fig.  15.  The 

data  show  that  there  are  considerable  qualitative  differences  between 

the  constant  density  and  variable  density  correlations.  It  could  be 

argued  that  proper  definitions  of  p*  and  p*  were  not  provided  by  Eqs. 

(43)  and  (44).  However,- a recent  study  (Ref.  25)  which  has  attempted 

to  correlate  constant  density  and  variable  density  combustor  flows  using 

the  C*  parameter  but  with  different  definitions  for  p*  and  p*  had 
i k o 

inconclusive  results.  Therefore,  at  present,  suitable  definitions  for 
pj^  and  p*  to  correlate  constant-density  and  variable-density  combustor 
flows  do  not  exist.  Hence,  additional  understanding  appears  necessary 
in  order  to  formulate  a useful  variable-density  confined  mixing  simili- 
tude parameter  such  as  C*. 


5.0  CONCLUSIONS 

An  analysis  of  the  literature  revealed  that  before  confident 
design  of  combustor  processes  based  on  recirculating  flows  can  be  pro- 
vided, much  more  experimental  and  theoretical  study  of  ducted,  axisym- 
metrlc,  turbulent  jet  mixing  with  recirculation  is  required.  Detailed 
recommendations  for  future  study  and  research  of  such  flows  are  pre- 
sented in  Ref.  7.  The  present  investigation  obtained  a set  of  detailed 
experimental  data  defining  ducted,  coaxial,  variable-density  turbulent 
jet  mixing  with  recirculation.  The  primary  or  central  jet  was  air,  and 
the  annular,  secondary  stream  was  hydrogen.  The  study  provided  the 
axial  velocity,  hydrogen  concentration,  and  wall  static  pressure  fields 
for  three  sets  of  recirculating  flows  for  the  hydrogen-to-alr  velocity 
ratios  (0.02,  0.036,  and  0.07).  For  these  experiments,  the  density 
ratio  of  the  hydrogen  to  the  air  was  0.069,  and  the  mixing  duct  to 
jet  radius  ratio  was  10. 
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The  conclusions  from  the  present  Investigation  are: 

1.  The  velocity  fields  obtained  are,  in  general,  similar 
to  constant-density  recirculating  flow  experiments  with 
the  exception  that,  at  the  velocity  ratios  of  0.036  and 
0.07,  "bumps"  appeared  in  the  centerline  axial  velocity 
decay  curves. 

2.  The  locations  of  the  points  of  zero  velocity  in  the 

flow  field  are  difficult  to  measure,  and  for  this  reason, 
the  characteristic  flow  stagnation  point  locations  on  the 
wall  could  only  be  tentatively  defined.  In  the  main, 
however,  axial  velocity,  hydrogen  concentration,  and 
wall  static  pressure  fields  are  sufficiently  well  defined 
that  the  data  can  serve  for  design  of  recirculating  flow 
combustors  or  for  the  evaluation  and  development  of 
predictive  techniques  for  such  flows. 

3.  The  theoretical  study  of  the  variable-density, 
recirculating  combustor  flows  showed  that  care  must  be 
taken  and  considerable  effort  expended  to  obtain  reasonable 
calculations  of  such  flows.  Problems  that  remain  to  be 
overcome  include  the  computation  of  the  static  pressure 
fields,  construction  of  nonuniform  mesh  configurations 
within  which  routine  calculations  can  be  made  with 
conservation  of  species  assured,  specification  of 
vorticity  at  a protruding  corner  from  which  flow 
separates,  such  as  at  a nozzle  lip,  and,  finally,  the 
development  of  more  accurate  models  for  effective, 
turbulent  or  eddy,  viscosity  for  variable-density 

flows.  In  general,  the  numerical  solution  procedure 
studied  herein  could  not  be  considered  a satisfactory 
predictive  technique  for  engineering  design  calculations 
associated  with  such  flows. 
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4.  A variable-density  correlation  parameter  (C*) , 
analogous  to  the  Craya-Curtet  parameter,  was 
defined.  The  wall  stagnation  point  locations 
were  correlated  with  C*  for  the  three  u cases 
investigated.  The  correlation  appeared  quantitatively 
different  than  that  for  constant-density  flows,  and 
hence  the  utility  of  such  a parameter  for  correlation 
of  ducted,  recirculating  flows  could  not  be 
determined. 
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Figure  1.  Schematic  of  turbulent  ducted  mixing  with  recirculation. 
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b.  Computational  mesh  indicating  the  definition  of  aspect  ratio 
Figure  2.  Configuration  of  the  computational  mesh  and  the  boundaries 
of  the  region  encompassing  the  recirculating  flows. 
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Eddy  viscosity  distribution  for  confined  recirculating 
flows  which  includes  potential  flow  regions. 
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Heao  Axial  Velocity  (ft/sec) 

Note:  Solid  Lines  Drawn  Tbrougb  Experlaeotal  Data 

Dashed  Lines  Indicate  Estijaated  Profile  Shapes 

Small  O Theory,  Given  by  Eqs.  (30)  and  (31) 

Small  • Theory,  Given  by  Eqs.  (32),  (33),  and  (34)  (see  Fig.  3) 

a.  u = 0.02 

Figure  8.  Comparison  of  experimental  and  theoretical  radial  distributions 
of  mean  axial  velocity. 
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Mean  Axial  Velocity,  ft/sec 

Note:  Solid  Lines  Drawn  Through  Experimental  Data 

Dashed  Lines  Indicate  Estimated  Profile  Shapes 

Small  O Theory,  Ug  Given  by  Eqs.  (30)  and  (31) 

Small  • Theory,  -j.  Given  by  Eqs.  (32),  (33),  and  (34)  (see  EiK 


b.  u = 0.036 
Figure  8.  Continued. 
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Mean  Axial  Velocity,  ft/sec 

Note:  Solid  Lines  Drawn  Through  Experimental  Data 

Dashed  Lines  Indicate  Estimated  Profile  Shapes 

Small  O Theory,  Uj.  Given  by  Eqa.  (30)  and  (31) 

Small  • Theory,  Ug  Given  by  Eqa.  (32),  (33),  and  (3A)  (see  Fig.  3) 

c.  Concluded 
Figure  8.  Continued. 
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Mean  Axial  Velocity  (ft/sec) 

Note:  Solid  Lines  Drawn  Through  Experimental  Data 

Dashed  Lines  Indicate  Estimated  Profile  Shapes 

Small  O Theory,  u-g  Given  by  Eqs . (30)  and  (31) 

Small  • Theory,  Ug  Given  by  Eqs.  (32),  (33),  and  (34)  (see  Fig.  3) 


c.  u = 0.07 
Figure  8.  Concluded 
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Hydrogen  Mass  Fraction 


Note:  Solid  Lines  Drawn  Through  Experimental  Data 

Dashed  Lines  Indicate  Estimated  Profile  Shapes 

Small  O Theory,  Given  by  Eqs.  (30)  and  (31) 

Small  • Theory,  Given  by  Eqs.  (32),  (33),  and  (34)  (see  Fig.  3) 


a.  u = 0.02 

Figure  9.  Comparison  of  experimental  and  theoretical  radial  distributions 
of  mean  hydrogen  mass  fraction. 
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b.  u = 0.036 
Figure  9.  Continued. 
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Note:  Solid  Lines  Drawn  Through  Experimental  Data 

Dashed  Lines  Indicate  Estimated  Profile  Shapes 

Small  O Theory,  Uj,  Given  by  Eqs.  (30)  and  (31) 

Small  • Theory,  Given  by  Ec- . (32),  (33),  and  (34)  (see  Fig.  3) 


c.  u = 0.07 
Figure  9.  Continued. 
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Figure  9.  Concluded. 
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of  nondimensionalized  mean  axial  velocity  on  the  centerline  of 
the  combustor  simulator. 


Viscosity  model  with  potential  core 
Figure  10.  Concluded. 
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of  hydrogen  mass  fraction  at  the  wall  and  centerline. 
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b.  u = 0.036 
Figure  11.  Continued. 
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Figure  15.  Axial  locations  of  forward  and  rear  stagnation  points  in 
recirculating  combustor  flows  for  constant  and  variable 
density  jet  mixing  as  functions  of  the  Craya  Curtet  number. 
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Table  1.  Coefficients  of  the  Conservation  Equation  (Eq.  (12)) 


Function,  (p 

^4> 

Ip 

P 

0 

1 

1 

pr^ 

2 

2 

0) 

r 

r 

C 

1 

\I 

1 

Sc^ 

Table  2.  Source  Terms  of  the  Conservation  Equations 


Function 


'I' 


-Ui 


u 


i£ 

9r 


3 . 3p, 


= 0 (in  the  present  study,  flows  are 
nonreactive) 


C 
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Table  3.  Uncertainties  in  Measured  and  Calculated  Parameters  (0) 


Parameter 

Name  and 

Symbol,  0 

Percision 
Index,  Sp  , 
Percent  <rf 

n 

Bias 

Limit,  B- . 
F'ercent  or 

Uncertainty 

Percent  of 

Characteristic 

Value  of  e for 
Establishing  the 
Uncertainty, 

Air  Supply  Total 
Pressure 

il.O 

± 3.0 

± 5.0 

250  ps:g 

Hydrogen  Supply 

Total  Pressure 

il.O 

X 3.0 

z 5.  0 

100  psig 

Air  Supply  Total 
Temperature 

iO.  15 

t 1.8 

t 2.  1 

530‘R 

Hydrogen  Supply- 
Total  Temperature 

zO.  IS 

i 

. 6 

± 2.  1 

SSO'R 

Primary  Nozzle 
i*lenum  Pressure 

tl.  0 

t 1 

.0 

t 5.0 

5 p«ig 

Secondary  Nozzle 
Plenum  Pressure 

i2.0 

i 6.0 

xlO.O 

1 ps.g 

Primary  Nozzle 

Throat  Static 

Pressure 

il.O 

± 3.0 

2 5.  0 

1 psig 

Kxhaust  Duct 

Static  Pressure 

t2.  0 

X 6.0 

ZlO.O 

1 psig 

V*  all  Static 

Pressure 

13.6 

X 3.6 

zl ..  4 

0.47  psig 

Probe  Pitot 

Pressure. 

H.gb  Range 

±1.0 

X 3.0 

2 5.0 

3 peig 

Probe  Pitot 

Pressure, 

L<iw  Range 

z2.0 

1 4.0 

2 6.0 

0.  025  ps:g 

Wall  -Static 

Pressure  for 

Velocity  Calcula- 
tions 

±2.  0 

2 4.0 

± 6.0 

0.  025  psig 

Hydrogen  Gas 
Concentrations, 

High  Range 

±4.0 

± 2.0 

210.0 

0.90 

Hydrogen  Gas 
Concentrations. 

Low  Range 

±1.0 

± 1.  5 

t 4.0 

0.  1 

Velocity,  u 

High  Range  at 

Cc  • 0.  1 

±0.  7 

± 

. 9 

1 3..1 

550  ft/sec 

Velocity,  u 

Low  Range  at 

C(^  0,90 

±6.9 

±14.  1 

±27.  9 

80  ft  sec 

Air  Supply 

Mass  Flow 

±1.6 

2 6.  7 

0. 07  Ibm  sec 

Hydrogen  Supply 

Mass  How 

il.6 

± ^-.7 

0.035  Ibm  sec 

Probe  Radial 

Position 

±0.  6 

t 6.2 

5 in. 

Axial  Distance 
between  T’robe 

T ip  and  Nozzle 
Assembly 

to.  2 



t 

.0 

1 1.4 

50  m. 
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Appendix  A contains  a summary  of  the  experimental  data  obtained  in 
the  present  study.  The  data  are  the  measurements  of  time-averaged  axial 
velocity,  hydrogen  gas  mass  fraction  in  the  interior  of  the  duct  and  on 
the  duct  wall,  and  the  wall  static  pressures  for  the  primary  and  secondary 
influx  conditions  shown  in  Table  4. 
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Table  A-2.  Mean  Hydrogen  Mass  Fraction  Data  - Interior  Field  Data 
a.  Mean  Hydrogen  Mass  Fraction,  Cc,  from  Test  A 


3.8  0.02 


0.007  1 0.021 


0,388  0.362 


0.191  0.216 


171.7  0.135 


0.  153  0.  152 


56.3  I 0.137 


5 

0.  14 

0 

0.  143 

0.035  0.390  0.516  0.515  0.408 


38.2  0.047  0,037  0.082 


0.207  0.206 


0.206 


0.20!)  0.210 


0. 66  1 

0.  77 

0.  85 

0.457 

0.  448 

0.  449 

0,  255 

0.  257 

0.  255 

0.  221 

0.219 

0.211 

0.211 

0.  212 

0.  207 

0.  206 

0.  208 

0.208 

IQQQI 

0.  207 

0,  209 

0.  209 

Table  A-3.  Hydrogen  Mass  Fraction  at  the  Wall 


Axial  Distance, 
x/R 

P 

Test  A 

Test  B 

Test  C 

Test  D 

Test  E 

Test  G 

0 

1.43 

0.  371 

0.334 

0.  683 

0.  624 

0.  772 

0.  986 

5.25 

0.  383 

0.  332 

0.  713 

0.  646 

1.  00 

... 

9.07 

0.  272 

... 

0.  581 

0.  552 

1.00 

0.  872 

12.88 

0.  272 

0.  269 

0.  502 

0.  518 

1.00 

0.  958 

16.  70 

0.  237 

0.  244 

0.  494 

... 

0.  700 

0.  907 

20.  52 

0,  207 

0.  201 

0.  417 

0.442 

0.  993 

0.  742 

24.33 

0.  186 

0.  187 

0.  383 

0.  389 

1.  00 

... 

28.  15 

0.  169 

0.  188 

0.  346 

0.  326 

0.  937 

0.  763 

31.  97 

0.  159 

0.  156 

0.  314 

0.  306 

0.  832 

0.  698 

35.  78 

0.  159 

0.  141 

0.273 

0.275 

0.  811 

0.  655 

39.  60 

0.  149 

0.  137 

0.  256 

0.  267 

0.  564 

0.  627 

43.42 

0.  143 

0.  145 

0.  239 

0.  245 

0.  569 

... 

4 7.23 

0.  135 

0,  126 

0.228 

0.225 

0.535 

0.496 

51.05 

0.  135 

0.  128 

0.  241 

0.  239 

0.  463 

0.  462 

54.87 

0.  143 

0.  115 

0.  220 

0.225 

0.  486 

0.  444 

58.  68 

0.  139 

0.  130 

0.217 

0.228 

0.  405 

0.  465 

62.  50 

0.  136 

0.  134 

0.  216 

0.  223 

0.  444 

... 

66.32 

0.  130 

0.  137 

0.  206 

0.  208 

0.  417 

0.  409 

70.  13 

0.  130 

0.  131 

0.  209 

0.207 

0.  373 

0.  399 

73.  95 

0.  146 

0.  132 

0.  204 

0.211 

0.415 

0.409 

77.  77 

0.  133 

0.  128 

0.  204 

0.  225 

0.323 

0.411 

81.58 

0.  132 

0.  133 

0.200 

0.211 

0.  374 

... 

85.  40 

0.  126 

0.  137 

0.  197 

0.  200 

0.357 

0.  363 

89.22 

0.  128 

0.  126 

0.202 

0.  207 

0.  340 

0.  358 

93.03 

0.  137 

0.  131 

0.  196 

0.  207 

0.  340 

0.  357 

96.  85 

0.  134 

0.  126 

0.  199 

0.215 

0.  335 

0.095 

100,67 

0.  143 

0.  141 

0.  215 

0.222 

0.  380 

... 

104. 49 

0.  127 

0.  140 

0.  196 

0.  199 

0.  344 

0.333 

108. 30 

0.  129 

0.  127 

0.  202 

0.204 

0.329 

0.  345 

112.  12 

0.  138 

0.  130 

0.  199 

0.203 

0.333 

0.  326 

115.  94 

0.  135 

0.  128 

0.  204 

0.  214 

0.  340 

0.366 

119.  75 

0.  135 

0.  137 

0,  203 

0.  206 

0.  364 

... 

120.  71 

0.  324 

121.  66 
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Table  A-4.  Wall  Static  Pressure  Distributions 


Axial  Distance, 
x/R 

P 

Ap,  psf 

u = 0.02 

u = 0.036 

u * 0.07 

0 

0 

1.  43 

0 

5.  25 

0 

9,07 

0 

12.  88 

0 

16.  70 

0 

20.  52 

0 

24.  33 

0 

28.  15 

-0.  55 

0 

0 

31.97 

-0.  55 

0 

-0.  53 

35.  78 

-1,00 

-0.57 

-0.53 

39.  60 

-1.00 

-0.  80 

-0.  53 

43.  42 

-0.  55 

-0.  80 

-1.05 

47,  23 

0 

-0.  80 

-1.05 

51,  05 

-0.  55 

-0.  80 

-1.05 

54.  87 

1.  55 

0 

-1.05 

58.  68 

3.  10 

0.  85 

0 

62.  50 

4.  15 

1.  841 

0 

66.32 

5.  70 

3.31 

1.05 

70.  13 

6.  75 

4.  34 

1.05 

73.  95 

6.  75 

4.  63 

2.  10 

77,  77 

7.  30 

5,  67 

2.  63 

81.  58 

7.  76 

6.  19 

3.  33 

85.40 

7.  76 

6.  99 

4.  16 

89.  22 

7.  76 

7.  32 

4.  64 

93.03 

7.  76 

7,  32 

5.21 

96.  85 

7.  76 

7.  56 

5.  69 

100.  67 

8.  26 

7.  74 

6.21 

104.  49 

8.  26 

7.  74 

6.  78 

108.  30 

7.  76 

7.  74 

6.  78 

112.  12 

7.  76 

8.03 

7.00 

115.  94 

7,  76 

7.  56 

6.  78 

120. 71 

8.  26 

8.03 

7,31 

123. 56 

7.  76 

7.  56 

7.00 

127,  39 

7.  76 

7.  56 

7,00 

131.20 

7.  76 

7.  56 

6.  74 
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NOMENCLATURE 

A Area 

Bias  of  the  parameter  ((p) , units  of  iji 
Hydrogen  gas  mass  fraction 
Craya-Curtet  number,  (Eq.  (2)) 

C*  Variable  density  Craya-Curtet  number  (Eq.  (39)) 

C Transported  scalar  quantity 

2 

D~  Effective  diffusion  coefficient,  ft  /sec 

^E 

d^  Source  term  (Eq.  (12)) 

L Length,  ft  or  in.  as  noted 

Thickness  of  primary  nozzle  lip  in  the  theoretical  solution 
procedure  (Fig.  2) 

N Reynolds  number,  N = puL/p 

R6  Rs 

N Effective  Schmidt  number,  N = p /pD^- 

bcg.  Sc  E Cg 

Pj.  Total  pressure,  psfa 

p Static  pressure,  psfa 


AEOCTR-76  152 


Ap 


R 

P 


R 

s 


R 


u 


u 

m 


u 


u 

P 


u 

s 


Wall  static  pressure  difference  defined  as  the  wall  static 
pressure  at  a given  axial  distance  from  the  primary  nozzle 
exit  plane  minus  the  wall  static  pressure  in  the  primary 
nozzle  exit  plane,  Ap  = p(x)  - p(o),  psfa 

Radius  to  nodal  point  in  the  theoretical  solution  procedure 
Primary  nozzle  radius,  ft  or  in.  as  noted 
Duct  radius,  ft  or  in.  as  noted 

Radius  of  primary  nozzle  outer  lip  in  the  theoretical  solution 
procedure  (Fig.  2) 

Radial  coordinate 

Nondimensional  radial  coordinate,  r^  = t/R^ 

Precision  index  of  the  parameter  ($),  units  of  $ 

Total  temperature,  °R 
Time-averaged  axial  velocity,  ft/sec 

Maximum  axial  velocity,  also,  axial  velocity  on  the  centerline 
of  the  duct,  ft/sec 


Velocity  ratio,  u = u /u 

P s 


o 


Velocity  of  jet  from  primary  nozzle,  ft/sec 


Secondary  stream  velocity  at  given  axial  station,  ft/sec 
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u Secondary  stream  velocity  at  primary  nozzle  exit  plane,  ft/sec 

o 

Uj^  Characteristic  velocity  for  confined  jet  mixing,  ft/sec  (Eq.  (3)) 

u*  Characteristic  velocity  for  confined  Jet  mixing,  ft/sec  (Eq.  (4)) 

u*  Characteristic  velocity  for  variable  density  confined  jet 

mixing,  ft/sec  (Eq.  (42)) 

u**  Characteristic  velocity  for  variable  density  confined  jet  mixing 

ft/sec  (Eq.  (43)) 


core 


FS 


Uncertainty  in  the  parameter  ((fi),  units  of  $ 

Time-averaged  radial  velocity,  ft/sec 

Length  of  velocity  potential  core,  12  in  present  theoretical 
study 

Axial  distance  between  the  primary  nozzle  exit  plane  and  the 
forward  stagnation  point  on  the  duct  wall  (Fig.  1) 

Axial  distance  between  the  primary  nozzle  exit  plane  and  the 
rear  stagnation  point  on  the  duct  wall  (Fig.  1) 

Nondimensional  axial  coordinate,  x/R^ 

Axial  coordinate 

Coefficient  in  the  definition  of  corner  vorticity  (Eq.  (25)) 
Dynamic  viscosity,  Ibm/ft-sec 
Effective  dynamic  viscosity,  Ibm/ft-sec 
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I 


3 

p Density,  Ibm/ft 

Pj^  Characteristic  density  for  variable  density  confined  jet 

mixing  (Eq.  (44)) 

P*  Characteristic  density  for  variable  density  confined  jet 

mixing,  (Eq.  (45)) 

Ip  Dummy  parameter,  scalar 

p Stream  function  (Eqs.  (8)  and  (9)) 

n Vorticity  (Eq.  (11)) 

(iJ  Vorticity  divided  by  radius,  fl/r 

SUBSCRIPTS 

CA  Primary  nozzle  lower  lip  corner  point  (Fig.  2) 

CB  Primary  nozzle  upper  lip  corner  point  (Fig.  2) 

D Diameter 

E Effective  or  "Eastern" 

LAM  Molecular  property 

N Nodal  or  "Northern" 

o Zero  axial  station 

p Primary  jet  property 
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1 


r 1 

F 
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S "Southern"  or  secondary  stream  property 

W Wall  or  "Western" 


